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\Q , The analysis of singular regions in the NUT solutions carried out 
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in the recent paper (Manko and Ruiz, 2005 Class. Quantum Grav. 



qh 22, 3555) is now extended to the Demiahski-Newman vacuum and 

electrovacuum spacetimes. We show that the effect which produces 
the NUT parameter in a more general situation remains essentially 
the same as in the purely NUT solutions: it introduces the semi- 
^ ' infinite singularities of infinite angular momenta and positive or neg- 

ative masses depending on the interrelations between the parameters; 
the presence of the electromagnetic field additionally endows the sin- 
gularities with electric and magnetic charges. The exact formulae 
describing the mass, charges and angular momentum distributions in 
the Demiahski-Newman solutions are obtained and concise general 
expressions P n = (m + i^)(ia) n , Q n = {q + \b){\a) n for the entire set 
of the respective Beig-Simon multipole moments are derived. These 
moments correspond to a unique choice of the integration constant 
in the expression of the metric function uj which is different from the 
original choice made by Demiahski and Newman. 

PACS numbers: 04.20.Jb 
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1 Introduction 



The Demmiahski-Newman (DN) stationary axisymmetric spacetimes are 
(i) the Kerr solution j2] endowed with the NUT parameter jSj (the pure 
vacuum case), and (ii) the Kerr-Newman solution [4: endowed with the NUT 
parameter and magnetic charge (the case of electrovacuum) . The canonical 
Papapetrou form of the general DN class was obtained in j^j as the simplest 
N — 1 specialization of the extended multi-soliton electrovacuum metric [H] . 

The presence of the non-vanishing NUT parameter makes the DN space- 
times asymptotically non-flat. This complicates the analysis of the respective 
mass and angular momentum distributions since, as was shown in [7] , the sin- 
gular semi-infinite regions developed by the NUT parameter have in general 
infinite angular momenta and non-zero masses. Moreover, the presence of 
the electromagnetic field requires additional consideration of the distributions 
of electric and magnetic charges in these singular regions. As an interesting 
aspect of the analysis of DN solutions one may consider establishing the 
correspondence between the multipole moments of the asymptotically NUT 
spacetimes calculated on a 3-manifold of trajectories of the timelike Killing 
vector with the aid of Simon's procedure |Hj and the physical characteristics 
of the entire four-dimensional metric. 

The plan of our paper is the following. In section 2 we consider the 
Ernst complex potentials and corresponding metric functions defining the DN 
spacetimes; here in particular we obtain concise formulae for the entire set 
of the Beig-Simon (BS) multupole moments associated with these solutions. 
Section 3 is devoted to the analysis of the mass and angular momentum 
distributions in the DN vacuum spacetimes, also known in the literature 
under the name Kerr-NUT metric. The more general case involving the 
NUT parameter and the electromagnetic field is analyzed in section 4. Lastly, 
concluding remarks are given in section 5. 
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2 The Ernst potentials, metric functions and 
BS multipole moments of the DN solutions 



It is commonly known that the stationary axisymmetric electrovacuum so- 
lutions of the Einstein-Maxwell equations are most easily presentable and 
analyzed if one uses, on the one hand, the Papapetrou canonical line element 

ds 2 = /~ VW + dz 2 ) + p 2 dip 2 } - f(dt - udip) 2 , (1) 

where the coefficients /, 7, u> are functions of the cylindrical coordinates p, 
z only, and, on the other hand, the Ernst formalism developed in the papers 
|Hl CHj according to which the knowledge of two complex potentials £ and $ 
satisfying the Ernst equations [TU] 

(Re£ + $$)A£ = (V£ + 2<SV$) V£, 

(Re£ + $$)A$ = (V£ + 2$V$)V$ (2) 

(a bar over a symbol denotes complex conjugation, V and A are the 3- 
dimensional gradient and Laplace operators, respectively) is sufficient for 
the reconstruction of the metric coefficients entering (|T|L together with the 
electric and magnetic components of the electromagnetic 4-potential. 

Demiahski and Newman originally derived their solutions by employ- 
ing some specific complex coordinate transformations not preserving Papa- 
petrou's form of the line element; the Ernst potentials defining the general 
class of DN metrics were obtained later in [H] using the extended multi- 
soliton electrovacuum solution 6J constructed with the aid of Sibgatullin's 
method ^TJ ^] ; these potentials have the form 

kx — m — i(ay + v) ^ q + ib 

kx + m — \{ay — v) ' kx + m — \(ay — v) ' 

1 / \ 1 / s 

x = — (r+ + r_), y=— (r+-r_), 
Zk Ik 

r± = J p 2 + (z ± k) 2 , k = \J m 2 + v 2 — a 2 — q 2 — b 2 , (3) 

where the arbitrary real parameters m, u, a, q, b can be associated, respec- 
tively, with the mass, gravomagnetic monopole (NUT parameter), angular 
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momentum per unit mass, electric and magnetic charges. In the absence of 
the electromagnetic field (g = & = 0<^$ = 0) one obtains from @ the 
expression for 8 defining the vacuum DN spacetimes. When the NUT pa- 
rameter v is equal to zero, one arrives at the black-hole solutions considered 
by Carter [T3] . 

The potentials (jSJ) can be used for the calculation of the corresponding 
BS multipole moments [UH^llE]- F° r this purpose it is advantageous to em- 
ploy the Hoenselaers-Perjes procedure [T£], recently rectified by Sotiriou and 
Apostolatos ^7j, which involves the axis expressions of the Ernst potentials. 
For p = 0, z > Re(hi) we have 

( \ -e( n ^ z-m-i(a + u) 



z + m — i(a — v) 



f(z) = $(p = 0,z) = g+ ,f -. (4) 

z + m — i{a — v) 

Then, passing to the potentials £ and 77 via the formulae 

and considering the expressions of the functions £ = z£ and rj = zrj in the 
limit z — > 00, we obtain 

(m + , . * (ia) n 



_ (g + i&)z , ^-n^( ia ) n (r\ 

Remarkably, it turns out that the coefficients in expansions © coincide 
exactly with the actual BS multipole moments because, as one trivially ver- 
ifies, all quantities My, Qy, <%, #y defined by formulae (23) of [T7] (they 
describe the deviations of the coefficients in (JHJ) from the multipole moments) 
are equal to zero identically. Hence we arrive at the following elegant explicit 
expressions for the complex multipole moments: 

P n = M n + U n = (m + iz/)(ia) n , Q n = E n + iB n = (g + i6)(ia) n , (7) 
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where the real quantities M n , J n , E n , B n are, respectively, the mass, angular 
momentum, electric and magnetic multipole moments of the DN solution 
(J3J). Formulae (J7|) generalize in a very natural and straightforward way the 
Sotiriou-Apostolatos result ^7] derived for the Kerr-Newman spacetime. 
The explicit expressions for M n , J n , E n and B n are readily obtainable from 
(0) by considering separately the even and odd moments: 

M 2k = (-l)W fc , M 2k+1 = (-l) fe+ W fc+ \ 
J 2k = (-l)W fc , J 2k+l = (-1) W fc+1 , 
E 2k = (-l) k qa 2 \ E 2k+1 = (-l) k+1 ba 2k+1 , 
B 2k = (-l) k ba 2k , B 2k+1 = (-l) k qa 2k+ \ k — 0,1,... (8) 

We find it appropriate to make now several comments on the multipole 
moments obtained. First of all, one should remember that the BS moments, 
like the Geroch-Hansen multipole moments [T%l ll9|. are only well-defined for 
the asymptotically flat spacetimes. In the stationary axisymmetric case these 
are spacetimes whose metric coefficients in have the following behavior 
at spatial infinity: 

/-l, 7-0, u;-0. (9) 

Although the last condition on u does not enter explicitly into the proce- 
dure for the calculation of multipole moments because the latter are defined 
on a specific 3-dimensional manifold, it is normally taken into account by 
demanding that the NUT parameter of the solution be equal to zero. As 
was already remarked by Simon [B] , the presence of the NUT parameter and 
magnetic charge may cause the three- and four-dimensional descriptions of 
a stationary electrovac solution topologically incompatible. 

Secondly, the multipole expansions define the corresponding Ernst po- 
tentials uniquely, which is a consequence of the mathematical theorems on 
analytic functions. On the other hand, the corresponding metric functions 
7 and uo are constructed from the potentials £ and $ up to two arbitrary 
real additive constants. The constant in 7 can be easily adjusted to have 
the desired behavior of 7 even in the presence of the NUT parameter; be- 
sides, 7 is involved explicitly in the 3-dimensional manifold which is used for 
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the calculation of the multipole moments and hence can be always defined 
uniquely by imposing the asymptotic flatness condition. The unique choice 
of the integration constant in the expression of uj leading to uj — > at spatial 
infinity can be realized only in the absence of the NUT parameter. When the 
parameter NUT is present in the solution, the formal multipole expansions 
of the corresponding Ernst potentials are still possible, but the respective 
4-dimensional spacetime is not globally asymptotically flat for any choice of 
the integration constant in uj. Therefore, an important non-trivial question 
arises of whether the multipole moments of the asymptotically NUT space- 
times obtainable from the Ernst potentials describe correctly the sources, and 
which is the precise choice of the integration constant in uj corresponding to 
those moments? 

Thirdly, because it is clear that the analysis of sources in the asymp- 
totically NUT solutions cannot be restricted to only the calculation of the 
multipole moments, one also needs a procedure of the evaluation of physical 
quantities which would be applicable to asymptotically non-flat spacetimes. 
The Komar integrals |20j then seem to be the best option as they had al- 
ready proved to be very efficient for treating the black-hole properties in the 
external gravitational fields j2U 122] • In our recent paper j2j we have applied 
the Komar integrals to the analysis of the mass and angular momentum dis- 
tributions in the NUT spacetimes; the main result of that paper consists 
in establishing the unique value of the integration constant entering the ex- 
pression of uj at which the total angular momentum is a finite quantity and 
hence is in accordance with the multipole structure. One may expect that an 
analogous criterion of the choice of the integration constants can be worked 
out on the basis of Komar integrals for the DN spacetimes too. 

Let us proceed with the description of the DN spacetimes in the Ernst 
picture by writing out the corresponding metric functions of the electric A± 
and magnetic A3 components of the electromagnetic 4-potential worked out 
in 0: 

= K 2 {x 2 -l)-a 2 {l-y 2 ) g27 = K y -l)-a 2 (l-y 2 ) 
(kx + m) 2 + [ay — u) 2 K 2 (x 2 — y 2 ) 
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uj = 2v(y + Ci) 



a(l — y 2 ) [2(rriKx — vay + m 2 + v 2 ) — q 2 — b 2 ] 



k 2 (x 2 — 1) — a 2 (l — y 2 ) 

q(nx + m) + — ay) 



^3 = &(C 2 - y) + 



(kx + m) 2 + (ay — v) 2 ' 
(1 — y)(ay + a — 2v)[q(Kx + m) + b(v — ay)] 



(10) 



(kx + m) 2 + (ay — v) 2 



Here Ci and C2 are two arbitrary real (integration) constants whose concrete 
values corresponding to the multipole moments (j2J) cannot yet be pointed out. 
Formulae fllOD can be formally called the entire family of DN spacetimes, and 
it is worthwhile mentioning that the particular choice of C\ and C2 made in 
jH] does not lead to the geometry defined by the multipoles (JJJ). 

It is convenient and instructive to consider the pure vacuum DN space- 
times separately from the electrovacuum ones (as was done in the original 
Demiahski-Newman paper 0): the simpler, vacuum case provides a solid 
basis for understanding the general structure of the DN sources, and we are 
passing to its consideration in the next section. 

3 The mass and angular momentum distri- 
butions in the DN vacuum subclass 

The vacuum subclass of DN spacetimes is normally known in the literature 
as the "genuine" Demiahski-Newman metric. It is also known under the 
name "(combined) Kerr-NUT metric" given to it by Demiahski and New- 
man themselves because it contains one additional (NUT) parameter com- 
pared to the Kerr solution 2j. Using formulae of the previous section let 
us write out the Ernst potential £ (the other potential $ is equal to zero 
identically) and corresponding metric functions determining these particular 
vacuum spacetimes: 



8 = f + in 



kx — 771 — i(ay + v) 



2(ukx + may) 




(kx + m) 2 + (ny — v) 2 
k 2 (x 2 -\)-a 2 (\-y 2 ) 
k 2 (x 2 — y 2 ) 



2 ■ 



f 
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2a(l — y 2 )(rriKX — vay + m 2 + i/ 2 ) 
a; = 2z/(y + C) 



K 2 (a; 2 — 1) — a 2 (l — y 



2^ 



x = ^-(r + + r_), y = — (r+-r_), r± = yV + (2 ± k) 2 , 

k = "\Aw 2 + z^ 2 — a 2 , (11) 

where the real constant C\ from (fTU|) is now called C, while the expression 
for k does not contain the parameters q and b (cf. formula Q of the previous 
section). 

When C / ±1, i/ / 0, the geometries (jllj) are characterized by two semi- 
infinite singular sources located on the symmetry z-axis: p = 0, z > Re(ft) 
(the upper singularity), and p = 0, z < — Re(/c) (the lower singularity). In 
the particular cases C = ±1, v ^ 0, only one semi-infinite singularity is 
present, which is analogous to the pure NUT case [Zj- If k is a real non-zero 
quantity, i.e. m 2 + v 2 > a 2 , the symmetry axis is divided into three parts: 
the upper region I (k < z < oo), the intermediate region II (\z\ < k), and 
the lower region III (— oo < z < — k), see figure 1(a). For a 2 > m 2 + u 2 , 
K becomes a pure imaginary quantity (k = i\/a 2 — m? — v 2 \ and only two 
different regions, I and III, will then be present on the z-axis (figure 1(b)), 
the cut between the points — k and +k representing a superextreme central 
object. 

Since our major interest lies in the semi-infinite sources which are orig- 
inated by the NUT parameter, in what follows we shall mainly concentrate 
our consideration on the subextreme case m 2 + v 2 > a 2 and we will pay 
less attention to the superextreme case which can be considered as a sort of 
degeneration of the sources structure (the region /J of the symmetry axis 
disappears) . 

In the subextreme case the calculation of Komar quantities can be suit- 
ably carried out with the aid of Tomimatsu's formulae obtained in (2H1 : 

M = --u[Q(z= z 2 ) - Sl(z = zi)], 
J = --tu(z 2 - z x ) - Iuj 2 [Q{z = z 2 ) - Q(z = Zl )\, (12) 
where M and J define, respectively, the mass and angular momentum of the 
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part of the symmetry axis [zi, z 2 ] on which u takes the constant value. 

The functions involved in formulae (J 12)) are uj and Q taken on the sym- 
metry axis, and they have the following form there: 

Region I. This part of the z-axis is determined by p = 0, z > k, or 
equivalently by x — z/k, y — 1. Then, substituting these values of x and y 
into ^2 we obtain 

2(vz + ma) 

w = 2i/(l + C7), n = - ? ,\ 2 , ( ' Vr 13) 

(z + m) z + (a — v) 1 

Region II. The intermediate region which is associated with the central 
rotating body is defined by p = 0, \z\ < k, or x = 1, y = z/k, thus providing 
us with the following values for u and Q: 



2 r , . 2n(maz + k 2 i/) 
- m (« + m + i/ aC + v , = v ^ 



Region III. Here p = 0, z < — ac, i.e. x = —z/k, y = — 1, and we have 

W = 2„(C-1), n= 2( 7 + 7 a > (15) 

(2; — m)^ + (a + z/)^ 

We first apply formulae (JT2J) to (fI3 )) -([T5 J) for getting the masses Mi, M 2 , 
M 3 of the corresponding regions I, II, III, the result is 

■,(l + Q(q + „) M 2 = m+^ C ' + !/) 



2(ac + m) k + m 

A 2{K + m) y 1 

It is easy to verify that the total mass M = Mi + M 2 + M3 is equal to m. 

From (|T2j) it follows that the entire angular momenta of the semi-infinite 
singular regions are infinitely large quantities since in the case of the region I 
one has to put z 2 = +00, and in the case of the region III z\ = —00. 
Therefore, as in the paper [Jj, we shall calculate the angular momentum 
Ji(zo) of the part k < z < z of the upper singular region, the total angular 
momentum J 2 of the central body (\z\ < k), and the part — z Q < z < —k of 
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the lower semi-infinite singularity. Formulae (|12|) then give 



Mzo ) 



v{\ + C) 
2 

J 2 = (a + Cv) 
1/(1 -C) 

^0 



i/(l + C)(a + i/) 



k + m 
v{aC + v) 

k + m 
y(l -C)(a- tj 
k + m 



+ 



i/ 2 (l + C) 2 (ra + ma) 
(z + m) 2 + (a - v) 2 



v 2 (l-C) 2 (vz -ma) 
(z + m) 2 + (a + v) 2 



• (17) 



The total angular momentum of the part \z\ < zq of the z-axis thus takes 
the form 1 

J(z = —Cu(z — 2m) + ma 

v 2 {\ + C) 2 (vz Q + ma) v 2 {\ - C) 2 (uz - ma) 



(18) 

(zo + m) 2 + (a — v) 2 (zo + m) 2 + (a + v) 2 
and one can see that the only choice of the constant C leading to the finite 
value of J in the limit zq — > 00 is C = 0. In particular, in the original DN 
solution characterized by C = —1 the only (lower) semi-infinite singularity 
carries an infinite angular momentum. 

In the general case, the aggregate mass of two semi-infinite singular 
sources is equal to 

u(aC + v) 



M, 



M1 + M3 



(19) 



0, M agg 



k + m 

and it can assume either positive or negative values. When C 
is a positive quantity for all non-zero m, v and a. However, even in the 
cases when M agg is positively defined, the negative mass is present in the 
singular semi-infinite sources. Let us illustrate this by taking as an example 
the original DN vacuum solution (C = — 1). In this case M agg = M3 = 
v{a—v) I (n+m), so that M agg is positive for instance when v > 0, v < a < m. 
In figure 2 we have plotted the mass M 3 (z ) of the part — z < z < —k of 
the semi-infinite singularity, namely 



M 3 {z ) 



v(a — v) 2v(vzq — ma) 

k + m (zq + m) 2 + (a + v) 2 



(20) 



1 Thcre is a misprint in the formula (15) of \7\: the first expression in parentheses on 
the right-hand side does not contain a and should be read as (zq — 2m). 
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as the function of zq, for the particular choice of the parameters m = 5, 
v = 3, a = 4. It follows that M 3 (z ) is a monotonously increasing function 
on the interval 4.243 < zq < 20.272, taking its maximal value 0.681 at 
Zq = 20.272, and it is a monotonously decreasing function on the interval 
20.272 < zq < oo, reaching asymptotically the value 0.325, which means 
that the latter interval is entirely composed of the negative mass. 

Moreover, even in the special case v = a when the total mass of the 
DN semi-infinite singularity is equal to zero, the qualitative picture of the 
mass distribution in the singularity is similar to that shown in figure 2: the 
interval with the positive mass is followed by the distribution of the negative 
mass. This is illustrated in figure 3 for the particular choice of the parameters 
m = 5, v = 3, a = 3, the total mass of the singularity reaching asymptotically 
the value M 3 = 0. 

When the central body is a superextreme object, i.e., in the case a 2 > 
m 2 + v 2 , the structure of singularities preserves the main characteristic fea- 
tures of the subextreme case: the aggregate angular momentum of the sin- 
gularities assumes a finite value only when the constant C is equal to zero 
(for the non-vanishing m, v and a), and the semi-infinite singular regions 
involve negative masses. In figure 4 we have plotted the mass distribution in 
the singularity of the original DN vacuum solution [C = —1) for two partic- 
ular parameter sets which have common values m — 2, v — 1 and differ in 
the value of a: in the first set a = 4 (figure 4(a)), and in the second a = —4 
(figure 4(b)). 

Since only the case C = is consistent with the multipole moments de- 
fined by the Ernst potential below we shall give the general expressions 
for the masses and angular momenta exclusively for this superextreme case, 
the subindex 2 referring to the central object: 

mav , , m(m 2 + v 2 + a 2 ) 2 

M 1 = -. M, 



M 



m 2 + (a — v) 2 ' {m 2 + v 2 + a 2 ) 2 — 4a 2 z/ 2 

mav 



3 



m? + (a + z/ 



2 ' 



vzq mav 2 v 2 (vzq + ma) 



Jl( z o) — o 2 i I V2 + 



m 2 + (a — v) 2 (zq + m) 2 + (a — v) 2 ' 
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2mais 2 (m 2 + v 2 + a 2 ) 

J 2 = ma + 



(m 2 + v 2 + a 2 ) 2 — 4a 2 



2 ' 



, / x _ vzq mav 2 v 2 {vzq - ma) 

j3[Zo) ~ 2 m 2 + (a + z/) 2 (z + m) 2 + (a + v) 2 ' [ ' 

As before, M\ and M3 are the total masses of the upper and lower singular 
regions, M2 and J2 are the mass and angular momentum of the central object, 
while J\{zq) and J^{zq) denote angular momenta of the parts < z < zq and 
—z < z < of the singular regions, respectively. Note that since formulae 
(12 lj) have been derived under the supposition a 2 > m 2 + z/ 2 , the limit a —>■ 
in (|2*Tj) is impossible, as m and v are real quantities by definition. 

It is easy to see from (|2*T]) that the aggregate mass of two singularities 
always takes a negative value for any positive m and v 7^ 0: 

Ama 2 iy 2 

Mi + M 3 = — (22) 
(m z + v A + a z ) 2 — 4a z u z 

In figure 5 we have plotted the mass distributions in the semi-infinite singu- 
larities for a particular choice of the parameters m — 1, v — 2, a — 3. 



4 The electrovacuum case 



In the presence of the electric and magnetic charges the general picture of 
the singular regions due to the NUT parameter is similar to that of the 
vacuum case. However, the singularities will now be electrically and mag- 
netically charged, and their masses and angular momenta will have both the 
gravitational and electromagnetic contributions. Since the properties of the 
semi-infinite singularities are practically independent of the central body, we 
shall restrict our consideration to only the subextreme case defined by the 
inequality m 2 + v 2 > a 2 + q 2 + b 2 . The symmetry axis then is divided into 
three regions, as in figure 1: k < z < 00 (region I), \z\ < k (region II) and 
— k < z < —00 (region III). The functions u, Q = ImS, $ and A3 which 
are needed for the calculation of the Komar quantities show the following 
behaviour in the regions I, II, III: 
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Region I. For p = 0, K<z<oowe obtain from (jHJ) and (|l(Jj): 

2(vz + ma) 



oj{V) = 2u(d + 1), 



(z + m) 2 + (a — i/) 



2 ; 



$(D = - T g + f y A 3(1) = -6(C 2 + 1), (23) 

z + m — i(a — v) 

where the subindex 1 in parentheses (and the subindices 2, 3 below) signi- 
fies that the respective function is calculated on the indicated part of the 
symmetry axis. 

Region II. When p = 0, \z\ < k, we have 



(k + m) 2 + a 2 + v 2 „ 2n(maz + k 2 u) 

c 

iq) 



U(2) = 2C x v H , £1(2) = --, ^— — Yi — , \2> 

a (az — v) A + k z (k + m) z 

k*l n — i/7 

$(2) = 



— kz/ + i/«(«; + m) 
q(K + m) + ub 

A 3(2) = —<^20 

a 

[b(az — kv) — nq(n + m)] [(k + m) 2 + (a — v) 2 } 
a[(az — kv) 2 + k 2 (k + m) 2 } 

Region III. For p = 0, n < z < oo one gets 

2(vz + ma) 



(24) 



W( 3 ) = 2l/(Ci - 1), fi( 3 ) 



(2; — m) 2 + (a + z/) 2 ' 

^ q + ib 

$ ( 3 ) = 7 T v 

2: — m — i(a + z/j 

A 3 3 - -o(<^2 - 1J H 7 77j— 7 — ; — 77j— • (25) 

(2 — mr + (a + z/)^ 

The calculation of the masses Mj, electric Qi and magnetic 5j charges, 
as well as angular momenta Jj of the regions I — III can be carried out with 
the aid of Tomimatsu's formulae derived in the paper 



M { = Mf + Mf = — j ' - 2a; w Im($ (i) $ Wi ,)] ( fe 

1 f u * 



2J dl 
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— -IjJU)\VLu\(z = Ui) - 


Qu)(z = di)], 


Q l = 


iu; (i )Im[$ (i) (2: = u t ) - 


-<f>(i)(z = di)], 


B t = 


--W(i)Re[$(i)(z = Ui 


= di)], 


Ji = 




(j)[2 + U(i)Q(i) jZ — 




1 f u i 

+ 2J d ^)^3lm($ (0 


,z)dz, 



(26) 

where the superscripts G and E denote the decomposition of masses and 
angular momenta into the gravitational and electromagnetic components in- 
troduced by Tomimatsu following Carter's paper |13j . 

The total masses of the regions I, II, III obtainable with the aid of the 
above formulae (jSEj) in which one has to set Ui = oo, d\ = w 2 = ^2 = ^3 = 
— k, ds = —00, have the form 

_ v{pi + l)(Kv + ma) 
(k + m) 2 + (a — v) 2 ' 
( 2uUu + ma) W 2av(C 1 - 1) \ 

y [K + my + (a — vy J \ [k + my + (a + / 
= Kgi ~ l)(KU-ma) . , 

3 (k + m) 2 + (a + v) 2 ' 1 J 

and it is easy to verify that the total mass M tot of the solution is simply 

3 

MM = Y / M i = m (28) 
i=i 

for all values of the constant C\. It is worthwhile pointing out that the 
cases C\ — — 1 and C\ — 1 correspond to vanishing of the upper and lower 
semi-infinite singularities, respectively. At the same time, for instance, the 
combination of the parameters kv + ma = 0, C% ^ — 1, which also leads to 
Mi = 0, does not annihilate the upper singularity, but only reflects the fact 
that the total gravitational and electromagnetic contributions into the mass 
of this singularity are equal in absolute values and have opposite signs; at 
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the same time, these contributions calculated for any finite interval of the 
singularity do not give zero. 

The distribution of the electric charge is described by the formulae 

v{d + 1)[6(k + m) + q(a - v)] 



Q 



(k + m) 2 + (a — v] 



2 



( 2u[b(K + m) + q(a - v)\ \ ( 2av(C 1 
Qi — \q 4 



(ac + m) 2 + (a — v) 2 j \ (k + m) 2 + (a + v 



\2 



v{Ci - 1)[6(k + m) - g(a + i/)] 



3 



i=l 



(k + m) 2 + ( r/ /') 

whereas for the distribution of the magnetic charge we obtain 
u(Ci + l)[q(K + m) - b(a - v)\ 



By 



(k + m) 2 + (a — v) 2 



B = / _ 2%(K + m)-6(a-^)] \ / + 2a^(d 



(k + m) 2 + (a — i/) 2 / \ (k + m) 2 + (a + u) 2 J 

B3 = _^-lMK + rn) + b (a + v)} ^ s ^ = 6 . (30 ) 
(k + m) 2 + (a + z/) 2 S 

One can see that the total charges too do not depend on the constant C\. 

The distribution of the angular momentum in the electrovac DN solution 
reminds qualitatively the distribution in the pure vacuum case, but the cor- 
responding formulae turn out to have a very cumbersome form and so will be 
not given here (the reader can work them out straightforwardly from (J2EJ))- 
We only comment that the choice C\ = is obligatory for having a finite 
value of the total angular momentum, and once this choice is made, the other 
constant Gi can be adjusted in such a way that the total angular momentum 
be equal to ma, in accordance with the multipole moments (J2J). 



5 Concluding remarks 

Like in the NUT solution, the semi-infinite singularities in the combined 
Kerr-NUT and electrovac DN spacetimes carry infinite angular momenta, 
and the only possibility for them to have a finite total angular momentum 
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is assigning zero value to the integration constant in the expression for the 
metric function u, in which case the two singularities will possess oppositely 
oriented angular momenta cancelling out the infinities in sum. The total 
mass of the singularities in this special case assumes a negative value. In 
contradistinction to the purely NUT solution where the counter-rotating 
singularities and a static central body form a system which is antisymmet- 
ric with respect to the equatorial plane, in the DN vacuum and electrovac 
spacetimes the non-zero parameter a makes impossible such an additional 
equatorial symmetry. The analysis carried out in the present paper and in 
[Zj shows how careful one should be when he tries to establish the multipole 
structure of spacetimes possessing the NUT parameter. At the same time, 
it is surprising how elegantly the additional parameters v and b of the elec- 
trovacuum DN spacetime generalize the known multipole expressions derived 
for the Kerr-Newman solution. 
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Figure 1: Different parts of the symmetry axis (a) in the case of the subex- 
treme DN spacetimes; (b) in the case of the superextreme DN spacetimes. 
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Figure 2: The mass distribution in the semi-infinite singularity of the subex- 
treme DN solution. The particular choice of the parameters is: m — 5, v — 3, 
a = 4, C = -1. 
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Figure 3: The mass distribution in the semi-infinite singularity of the subex- 
treme DN solution, the total mass of the singularity being equal to zero. 
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Figure 4: Dependence of the mass distribution in the original DN semi- 
infinite singularity on the angular momentum of the central body. 
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Figure 5: The mass distributions in the upper and lower singularities of the 
superextreme DN solution possessing finite total angular momentum (figures 
(a) and (b) respectively). 
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